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Using Momentum Transformation 
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We introduce a novel approach to model spatial discontinuities, such as thin surfaces and metasurfaces, in the 
momentum space (spatial spectrum). The approach is called "momentum transformation" to emphasize its spatial 


spectral nature and avoid confusion with temporal spectral techniques. 


The momentum transformation approach is analogous to expressing Maxwell's equations in the sense of distributions* 
but in the momentum space and is applicable to scalar and vector waves. The momentum transformation is particularly 


suitable for the scalar case due to its simplicity. 


The momentum transformation has been applied to metasurface synthesis, where it yields a complete description of the 
metasurface in the momentum space and provides physical insight into the transformation operation. 


These slides are adopted from an oral presentation at URSI GASS 2014, Beijing, China. 
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Conservation, Comparison with Optics, Momentum Transformation. 
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3D Metamaterials Metasurfaces Frequency Selective Surfaces 
(FSS) 






Soukoulis et al. 
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@ Signal Processing @ Reduced Dimensionality @ Filtering & Polarization 
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@ Bulkiness @ Extensive Processing @ Limited Signal Processing 
@ High Loss @ Low Loss @ Uniformity Restriction 


@ Fabrication Difficulty ®© Fabrication Easiness 
@ Incidence Flexibility 


Dimensional Reduction of Metamaterials & Functional Extension of FSSs , 
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Metasurface : 





Output Wave V 


Goals: 


e Systematic/universal 

e Scalar, vector, tensor 

* Infinite & truncated MS 
e PW & any other wave 

e Manipulation at distance 


e Easy to Implement 


enefits: 


SYNTHESIS — TWO-STEPS: 


* Determination of metasurface transfer function | ° Physical insight 
* implementation (scatterring particles) * Inspiring new devices 
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MT*: Momentum Transformation 
(К, not "spectral" avoid confusion with c) 
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тоо Foundation for the MT Method: Momentum Conservation (2. 


AMES 
M. A. Salem and C. Caloz, “Manipulating Metasurface 
light at distance by a metasurface using 
momentum transformation,” Opt. | 
Express, vol. 22, no. 12, pp. 14530- 
14543, Jun. 2014. 








| 
Input | Output 
f(z,y;0 ) JH | M seio? 
F. 2 
T UU 
Momentum qL L 
Transformation 


n(x, y) h(r.yimy)zh(r-z.y-y) 
со со Green function 
Linear shift-variance: 9 (2, y) - | J fe) hi үтен ау 
Zero EM thickness — locality: ^ (z', y';x,y) — ó(x' —x,y' — y)v(x',y) 
Consequence of locality: g (x, y) = f (x,y) n (x,y) 
Fourier transform: 9 (kz, ky) — f (kz, ky) EX hes, ky) 


Momentum shift-invariance - 2-0 ~ 
= momentum conservation: f (ka — Ку, Ку © Ку) g (k; — Ах, ky © Ку) 
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Locality: g (x, y) = n (x,y) f (£, y) > n (£, y) = m i 





1 
Ме{азиНасе: 7 (т, у) = 9(2,у)5(т,у), S(x,y) = 1 (т, у) 


Momentum conservation: f (kz, ky) = 9 (kx, ky) * € (kx, ky) 





Plane wave (Fourier) expansion at z = 0 of = /,”7,9: 


1 Too Too | 
Y (x,y) 2 — w (Ke, ky) etthett+ku) dk.dk Momentum wave equation: 

27 7 е 

a pes ke = ko + ko + ke 
(hee, ky) f° V (a, y) e Erk) dady 
NM (e t) 
UR. —ikzd „+ilkaz+kyy) 
Specification at z = d: g (2,7; 2 = d) = > | J Cha kyo Е а, 
Was 
ga (k; Ñy) 


Reverse _ _ _ 
propagation: g(k,, ky) > Jalka, ky) = 
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ие Momentum Conservation Method: Vector Case ес 





VPW orthogonality: J J Е (x,y; kz) x не (2, у; к.) атау - Py’) (5, - к.) 


k, eigenmode expansion: 


— 


Е (2, у; 2) = C* (k,)E V. V (x, y; k;) e'^*? dk, 


Lo, y, 0 
M" | М x H*VPW+T(r y; k,)dxdy 


& VPW orthogonality 





deu xVPW : 
С = Бур J. J = ) x HYYPW (x, y; kz) dedy 


expansion 


coefficients: (+ (kz) — mm] j Et (x,y) x H*VPW (z, y; k;) dzdy 
k — 00 J —oo 


Reverse propagator: 9 (kz, k,;k,d) = oik. d 


Scalarized (correspondence with scalar case) solution: ( (К„) ж й(К„) = CT (k.)@(k; d) 


Modal transformation, using k? = k? + ke + к2: 
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д. W. Goodman, Introduction to Fourier Optics, 2" edition, 1996. 
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System input system output 





f (2,7, Zin) CAT: g (£, Y, Zout) 
__ g (T, Y, Zout) 
t(x, y) = #2? Zout) 
P (2, y, Zin) 


EXAMPLE 4.1-1. Scanning. ^A thin transparency with complex amplitude transmittance 

I-A f(z,y) = exp(jraz*/Af) introduces a phase shift 2zó(r, y) where ó(z,y) — —2z^?/2Af, so that 

Not limited to the wave is deflected at the position (x,y) by the angles 6, - віп (Адф/дг) - віп ! (—z/f) 

апа 0, = 0. If |r/f| < 1, 0, = —r/f and the deflection angle 0, is directly proportional to 

the transverse distance z. This transparency may be used to deflect a narrow beam of light. If the 

transparency is moved at a uniform speed, the beam ts deflected by a linearly increasing angle as 
illustrated in Fig. 4.1-6. 


bulky elements 


Spatial Light 
Modulators 
(SLM) 





B. E. A. Saleh and M. C. Teich, Figure 4.1-6 Using a frequency- Figure 4.1-7 А transparency with transmit- 
Fundamentals of Photonics. 2007 modulated transparency to scan an optical tance f(z,y) — exp(jr1^/Af) bends the wave 
i 8 beam. at position т by an angle 0, 7: —z/ f so that it acts 


as a cylindrical lens with focal length f. 
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Sus a Transformation in MT and CAT бл 
Metasurface 
2 = (0) 

Input s 
J(z,y 0 )h, g(x, 430") 


RULIN 


Momentum transformation (MT): Complex amplitude transmittance (CAT): 


_ 9(%,у,0") оу 1902007) 
MOD = FE yO) OO Fy) 


Definition of input: 


1 (2, 9,0) =Л (т,у,07) +Л (т, у,0°) 


©) | . Роу 0) 
© ое Ішейгайіме Example: Anomalous Reflection & Refraction сага 


И P 
"А; 
4 (a, y, 07) _ 7 c ie sin(0* )a 
| 1 
ў (т, Y, 07) = cic sin(8*)z Metasurface 


z = Ü 
(passive, lossless) 


e.g. 2D problem, TE, wave 


Ө) | 
САТ Approach 
9 (2,5,0) 
t (x,y) = ——— — 
f(x,y, 07) 
7 eie sin(0*)z 
- cik sin(0!)x 


— r exp(ikz [sin(0*) — sin(0')]) 


Boundary conditions 
(2D-TE) 


f(z,y,07) 2 g(z, y, 07) 
д. f(x. y, z)l,g- = O; д(®, у, 2)|„—о+ 


f(x,y,z) = Ғіа,у,2)--Г(т,у,2) 


CAT & Boundary Conditions 


TE 
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What about the reflected wave? 


ik|sin(0*)x—cos(0*)z] 
А. 


7 c iF[sin(6* )z4-cos(0*) z] 
ošk[sin(0')z+cos(6')z] 
_ ešk[sin(0')z—cos(6'`)z] 


Cannot be used for synthesis: 


Not the required behavior 
Violating passivity 





р (т, у) = 9(т, у) — f'(x, y) 
тет sin(0°)z otk sin(0')x 


f(a, y. z) — 7 pik[sin(")x—cos(6")z] 


Е eik[sin(')x—cos(6")z 
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CAT Approach ik[sin(0! ) y —cos(0!) z] 


€ 
g (z,y,0*) 
t — _ _ — n 4, sin t 2; COS t Ж 
(ж, у) F(a, y,07) _reik[sin(6")«-+cos(6") A 


тете sin(0*)a Mum 
— "ziksin(8i)s- — ešk[sin(0* )z-r-cos(0*)z] 


= T exp(ikz [sin(0*) — sin(6')]) oik|sin(6')x+cos(6')z] 


Cannot be used for synthesis: 


Energy Conservation 
(same medium) 


f'(zx,y,07) 2 f'(x,y,07) 4 g(z, y, 07) 
— f'(z,y,0 ) t t(z,y)f (zy. 07) 


f'(,9,07) 2 (1— t(z,y)) f (z, 9,07) 


Not the required behavior 
Violating causality 





e - 
P EMQUE MT & Boundary Conditions 


MONTREAL 





MT Approach 


— g(2,y,07) 
ney) = f (x,y, 07) 


ro" sin(6t)x T a a. =. 
7 cik sin(0!)z Е ^yeik sin(0" )x 7 c ik[sin(0*)z--cos(0*) 2] 
eik[sin(6')x-+cos(6")z] 


Boundary conditions 
(in the sense of distributions) 


ныл 
0, g(z,y,z)|, 9 — 9s f(x, yz), g- = А (х, у) Correct physical description 


Related їо 77 (2, 7) 
А. (ж, у) = (7 (2,0) - 1) f(z; y, 0) 


f'(z,y,0 ) = g(x,y, 07) — f'(x,y, 07) — A(z, y) 
= туе?" sin(0* )x 


(©) 


m iti Poly >) 
245. ААТ & Boundary Conditions, Apparent Paradox W, 


4 L2 
RAMES 








Boundary conditions 
(in the sense of distributions) 


g(z,y,0*) — f(x,y,07) = " ” 


Momentum conservation (physics) 


j(k,, ky; 0*) f (kx, ky30~) = v) =m 5 (kx, ky;0T) — f (kz, ky;07) = 0 





Pagador? 


Discontinuous momenta 
(require distribution notation) 


T 


Particular solution: = Ü 
() 


Homogeneous solution: Supp À — су А = Л(х, у)0 (Ку, Ки) 


From momentum transformation: А = (m - 1) Í 


M. A. Salem and C. Caloz, Opt. Express, 22, 2014. 


Straigthforward generalization to vector case ! 


(Surface response function) 
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Output Wave 
(vortex hypergeometric Gaussian) з. 
| + 1ep/2 |ті ‚9 I | g(x, у) 
ue P M | C 
P (1-4 т] + 5) | р p ) | 
x —— i F [| -—=, т 1 — 
"(тау 2 ЕЯ 







`2 


Input Wave 
(ordinary Gaussian) 


Tic yx s әуе 
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POLYTECHNIQUE Delay-Start Airy Beam Generation Lion 






Output Wave (Airy beam with start at z = d) 


Al s - Sep (i13$ € 
Cu 56 С 
4/ 8072779 


(kz; k, d) = exp р 


8 = 7/20,6 -2/Ел 
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Input Wave 
(normally incident plane wave) 


f(z) = е 


g (z, z) = 
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O11 519 
n(x, y) => 
221 222 |(aeembedded) 


n(x, y;0) 


Rotation 


Scatterer topology 





. Rotation 


у/2. 








m(a, y; 0)| 


I— —I— + — + — 











0.4 


Rotatjon 


0.3 


0.2 


0.1 
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L] Metasurfaces: 
unprecedented EM control 





Q No universal synthesis 
— until now 
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nidirectional component 


_) Synthesis based оп 
transformation in momentum 
space for physical insight 


L1 Implementation using 
lookup maps 








L] Potential for myriads K. Achouri, M. A. Salem and C. Caloz, 
of applications "General Metasurface Synthesis Based on 
Susceptibility Tensors," arXiv:1408.0273. 
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